PUPT-2426 



Analytic quantum critical points from holography 

Jie Ren 

Department of Physics, Princeton University, Princeton, NJ 08544 
(Dated: November 21, 2012) 

We find exact, analytic solutions of the Klein-Gordon equation for a scalar field in the background 
of the extremal Reissner-Nordstrom-AdSs black hole. The Green's function near a quantum critical 
point for a strongly coupled system can be extracted holographically from an exact solution for the 
scalar at zero frequency (a;), but arbitrary momentum (k), mass, and charge. By examining the 
Green's function near lj = 0, there are two types of instability: the first one is triggered by a zero 
mode, and gives a hybridized critical point; the second one is triggered by the instability of the IR 
geometry, and gives a bifurcating critical point. The two types of instability can happen at the same 
time, and give a mixed critical point. Without tuning an extra parameter, only the second type of 
instability can happen at k = 0. At the critical point with the superfluid velocity, the scalar can 
develop either type of instability, depending on the parameters. The zero mode can also be obtained 
by tuning a double trace deformation. The phase diagrams can be analytically drawn. 
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I. INTRODUCTION 

Phase transitions for some strongly interacting sys- 
tems can be modeled by the gauge/gravity duality [I]- 
[3]. For example, holographic superconductors have been 
constructed in terms of asymptotic anti-de Sitter (AdS) 
spacetimes [11 [5]. In this paper, we consider quantum 
phase transitions, i.e., the phase transitions that happen 
at zero temperature. Understanding the quantum criti- 
cal points is a significant challenge in condensed matter 
physics, such as non-conventional superconductors. In 
the gravity description, a phase transition happens when 
a scalar field in an AdS background develops an instabil- 
ity. Solving the Klein-Gordon equation for the scalar field 
in the bulk gives the Green's function of a scalar operator 
in the boundary. The critical point can be identified by 
the non-analyticity of the Green's function at a; = 0. 

Previous studies have made major conceptual pro- 
gresses in understanding the quantum critical points in 
terms of the Reissner-Nordstrom (RN) black hole in AdS4 
[5] (for another system, see Ref. [7] ; the low-energy effec- 
tive field theory is also identified in Ref. f8|). The ex- 
tremal RN-AdS black hole has an AdS2 factor in its near 
horizon (hereafter IR for infrared) geometry This 
AdS2 factor plays an essential role in the properties of 
the system, and defines a universal intermediate energy 
phase [ini [H] . The behavior of the system near the quan- 
tum critical point is encoded in the Green's function near 
cj = 0, which contains both UV and IR data. The UV 
data can be solved from the Klein-Gordon equation at 
w = in the bulk, and the IR data can be analytically 
solved from the Klein- Gordon equation at arbitrary uj in 
AdS2 f3] . The zero mode, which triggers the onset of 
the instability of the system, belongs to the UV data, 
and relies on numerical calculations in AdS4. 

We find that if we use AdSs instead, the Klein-Gordon 
equation at w = can be analytically solved; the Green's 
function captures essential features of the RN-AdS sys- 
tem. We consider the standard/alternative quantization 



first. For a scalar field with mass m and charge q, the 
zero modes are solved as 

-. = ^-«±-^, (1) 

where n± is a nonnegative integer, A± is the scaling di- 
mension of the boundary operator, and Vk is the IR scal- 
ing exponent: 

A± 2 ± y/m^ + 4, (2) 

i^k := ^V^' + fc'- V+3. (3) 
2v3 

The zero modes are always at nonzero k, and the n± = 
mode triggers the onset of the instability. The instabil- 
ities and the corresponding quantum critical points are 
classified as follows (the name and interpretation of the 
quantum critical points are from Ref. [5]): 

• The first type of instability is triggered by a zero 
mode, which exists only if q is large enough. This 
instability gives a hybridized critical point, which 
is described by an order parameter in a Ginzburg- 
Landau sector hybridized with a strongly coupled 
sector, the CFTi dual to the IR AdS2. 

• The second type of instability happens when the 
IR scaling exponent Vk becomes imaginary, which 
implies the instability of the IR geometry. This in- 
stability gives a bifurcating critical point, for which 
the Green's function bifurcates into the complex 
plane. 

• The two types of instability can happen at the 
same time, and give a mixed critical point, such 
as a marginal critical point, which is described by 
a marginal term. 

We then study the quantum critical points at a; = 
and k = 0. Without introducing an extra parameter, we 
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can only have a bifurcating critical point. We can tune 
the superfluid velocity to reach all three quantum critical 
points. The phase diagram can be analytically obtained 
for the onset of the instability from the normal phase. 

Instead of the superfluid velocity, the zero mode can 
also be obtained by tuning another parameter k+ , the co- 
efficient of a double trace deformation in the CFT, giving 
a hybridized critical point 6, 12J. The analytic result al- 
lows us to draw the phase diagram for arbitrarily large 
TO^, where the numerical result is difficult to achieve. 
The Green's function (or susceptibility) for various crit- 
ical points can be obtained. 

This paper is organized as follows. In Sec. |llj we give 
the general solution of the Klein-Gordon equation at a; = 
0, and the Green's function near uj = 0. In Sec. |III[ 
we classify the instabilities according to the parameters 
and q, and give the critical values of the superfluid 
velocity. In Sec. |IV[ we consider the parameter of the 
double trace deformation, and draw the phase diagram. 
In Sec.jVj we conclude with some discussions. 



II. SOLUTION FOR THE KLEIN-GORDON 
EQUATION 

The extremal RN-AdSs black hole in Poincare coordi- 
nates i^ 



V /(2) 



/ = 1 - 



2z' 
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At = V6{l-z^), 



(4) 
(5) 



where x = {xi, X2, x^), and the gauge potential is A = 
A-tdt. We set the AdS radius L — 1, and the horizon is 
at z/i = 1. 

To obtain the Green's function for a scalar operator in 
the dual CFT, we will solve the Klein-Gordon equation 
for a scalar field <&. After the Fourier transform 



the equation of motion for 4> is 



(6) 



/ z 



(cj + qAt)'^ P 



P 



f 



= 0, 
(7) 

where we assume k = (fc, 0, 0) without loss of generality, 
q > 0, and is above the Breitenlohner-Freedman (BF) 
bound [13]: > mgp = —4. The horizon z = 1 is 
an irregular singularity for this equation. The in-falling 
boundary condition near the horizon z = 1 is 



W 



( ^) 

V Q{l-z)) 



(8) 



1 The action is S = J d^x^^{R- 12 - jF'^). 



z « 1 
AdSa X M3 
IR 



AdSs 
UV 



FIG. 1. The inner (near horizon) and outer regions, where 
the solutions of the Klein-Gordon equation are denoted by (pi 
and (j}o, respectively. 



where Wx^i^{z) is a Whittaker function. The asymptotic 
behavior near the AdS boundary 



^Az^-{l + ---) + Bz'^^l + ■■■). 



The retarded Green's function i^ 



(9) 
(10) 



When —4 < < —3, there is an alternative quantiza- 
tion, by which the Green's function is G = A/B [2]. 

To study the instability near a quantum critical point, 
we need to solve the Green's function near cj = 0. 
When it is sufRciently close to the extremal horizon, uj- 
dependent terms cannot be treated as small perturba- 
tions no matter how small uj is. In Ref. [3], a systematic 
method is developed for treating the extremal black hole 
system. We divide the geometry into inner and outer re- 
gions, as shown in Fig. [T] The inner region refers to the 
IR (near horizon) geometry, in which the Klein-Gordon 
equation can be exactly solved as Eq. ([s]). The outer re- 
gion refers to the remaining geometry, in which we can 
make perturbations for small uj. Then we need to match 
the inner and outer regions. 

In the inner region, the IR Green's function Gki'^) is 
solved as Eq. ( B6 1 (see Appendix [b]) 9] . In the outer 
region, the solution at small uj can be written as 



where 



(z) = ?7+(z) + gkiuj)ri_{z), 



(11) 



(12) 



At the leading order, the asymptotic behavior near the 
horizon z = 1 is 



T^f ^ [12(l-z)]-l/2±..^ 



(13) 



which we also use to fix the normalization of the solution. 
The asymptotic behavior near the AdS boundary z = 
is 



± 



(14) 



^ When A-t- — A_ = 2n, where n = 1, 2, ■ ■ ■ , there are extra terms 
bz^+ Inz (1 H ). When A+ = A_, ip = Az^\nz + Bz'^ H . 

^ We use the same normalization as in Ref. [6]. By another nor- 
maUzation, the Green's function is G = (2A_)_ — A)B/A. 
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The Green's function to the first order in is [S] 



G{uj,k) 



(0) , 



(15) 



Note that if we consider a neutral scalar, the first order 
terms in lo are zero, and then we need to expand to the 
second order. The analytic solution of ^ at w = be- 
low gives the leading order of the Green's function. By 
perturbation around = 0, we can obtain the higher- 



order coefficients. The Green's function can be general- 
ized to nonzero temperature when T « fi (chemical po- 
tential) by replacing the IR Green's function Gk{^) with 
Eq. (B7) in Appendix |b] The results we just described 
can be found in Rcf. [9 . We now go beyond Ref . [5] and 
obtain analytic solutions for a^^^ and b^^K 

When io — 0, we can solve </> in terms of hypergeometric 
equations. The general solution of </> for Eq. (FtI) at w = 
i|] 



(z) =Ciz 



(1 _ ^2)-l/2+., 



+C2Z 



(2z2 + l)-l/2+!.fc-HA_/2 

(1 - z2)-l/2+.. 
(2z2 + l)-l/2+!^fc+A+/2 



2F1 
2F1 



A_ - 1 
2 

A, -1 



q A_ - 1 
V3' 2 



3z' 



1) 



_q_ 

V3' 



vk -t- — ; A+ - 1; 



2z2 + l 



(16) 



The asymptotic behavior near the horizon z 1 is 



(1-z) 



{l~z) 



-l/2+Vk 



Ci 3-'^-/2r(A_ - i)r(-2i.fc) 



C2 3-^+/2r(A+-i)r(-2z.fe) 



' A_-l 



2 i^fe + ^'?)r(^^-z^fc 



^1 



A+-1 



C2 3-^+/2r(A+-i)r(2i/fe) 



2_)rf^ 



A, -1 



Ci 3-'^-/2r(A_ - i)r(2 



A_-l 



.r(^ + .. + -^)r( , 



73^ 



Vk 



73^ 



The asymptotic behavior near the boundary z — is 



^ Ciz'^- + C2Z^+ 



By Eqs. (|13| and (|14|), the solutions of of.'' and are 



(17) 
(18) 



/_18i/2+''^ • 3-^+/2p(^^ _ l)T{2vk) -l8i/2-^fc . 3-^+/2r(A+ - l)T{-2vk) \ 



\Jrr? + 4 



A, -1 



■A, -1 



181/2+^. . 3-A-/2r(A_ - l)r(2i/fc) 181/2-^. . 3-A_/2p(A_ _ l)r(-2j/fc) 



73)/ 



(19) 



It can be checked that a^^h^^ — a'^'fe'i^-' 
is satisfied. 



Vk/Vr. 



III. ANALYTIC GREEN'S FUNCTIONS 



A. Zero modes and the phase diagram 



* When A-^ is an integer, the two hypergeometric functions in 
Eq. are linearly dependent. We can choose another two 

linearly independent solutions as Eq. (lASl in Appendix rAl 



The zero mode is a gapless mode in the Green's 
function at w = 0. The singularity of a fermionic 
Green's function at a; indicates a Fermi surface, 
while the singularity of a bosonic Green's function at 
w = indicates instability. In the fermionic case, solv- 
ing G~^{ijj = 0, k) = gives k = kj?, where kj? is the 
Fermi momentum. Similarly, in the bosonic case, solving 
G^^{uj = 0,k) =0 gives the k = kg, where we use the 
subscript S for superfluid. The analytic solution enables 
us to solve for the zero modes, which are the only normal 
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modes. 

The zero modes are determined by a^'' = for the 

standard quantization, and b^^^ — for the alternative 
quantization. By Eq. (19 1, we have 



1 



V3 



(20) 



where n± is a nonnegative integer, the "+" sign is for 
the standard quantization, and the "— " sign is for the 
alternative quantization. Although the above result is 
derived under the assumption that A_|_ is not an inte- 



ger, Eq. (20) is generally valid when m~ > —4 for the 
standard quantization, and when —4 < < —3 for the 
alternative quantization (see Appendix [A| . The solution 
of k is denoted hy kg. The existence of a real kg requires 
> 0, which implies 



V3 



>n+ + 



A+-1 



(21) 



The first type of instability is initiated by the nodeless 
zero mode, n± — 0. We set n± = in the following. 
Solving k from Eq. ( 20 ) gives 



fc| = 12( 



q 



1 ± +4\2 



2q^ - 



q - — (1± +4) +2(1± v/m2~+4), 

3 



where 



g>^(i±y^ 



4). 



(22) 



(23) 



For the standard quantization, ks is always nonzero. For 
the alternative quantization, fcg is nonzero except for a 
special case: rr? = —3 and q = 0, in which, however, 
Eq. ([16]) is not valid and the solution of (j) by Eq. (44) 
below shows that k = is not a zero mode. Moreover, 
^5 is always nonzero for all n± > 0. 

If Eq. (23 1 is not satisfied, i.e., the zero mode does 



not exist, there can still be the second type of instability 
when z/fc is imaginary. The critical value of k is 



where 



fcj^R = 2q^ -771^-3, 



q^ > 



(24) 



(25) 



As explained in Ref. [9^, this instability is due to the 
backreaction of pair productions in the IR geometry for 
a charged scalar, or the mass below the AdS2 BF bound 
for a neutral scalar; the parameter set for which Vk is 
imaginary is call the oscillatory region. The imaginary 
IR scaling dimension 6^ = 1/2 + j/j. implies the confor- 
mality lost after the annihilation of two fixed points of 



the CFT dual to the AdS2, leading to an instability of 
the IR geometry [15l [16] . 

When a zero mode exists, there is always the second 
type of instability for k < km, by comparing Eqs. (23) 
and (1251): 



- < ^(1± v/m2+4). 



(26) 



and fciR < ks, by comparing Eqs. (22) and (24). The 



solution of fc as a function of and q is illustrated in 
Figs. J2I and [3] The blue lines are solutions of ks from 
Eq. (pOl); the outermost one is n = 0. The boundary of 



the shaded region is the solution of /cir from Eq. ( 24 ) . 

By perturbation around the solution at w = 0, we can 
obtain the retarded Green's function near the zero mode 



G = 



hi 



-U) — h2e'^'"'suj 



(27) 



where fcj^ — k ~ ks, and the quantities vs, hi, and h2 
can be calculated by the formulas given in Appendix C 
of Rcf . 12] . As explained in Ref. [3] , when the momentum 
changes from k > ks to k < ks, a. pole of the Green's 
function moves across the origin to the upper half com- 
plex oj-plane, signaling an instability. As a comparison, 
analytic fermionic Green's functions can be obtained for a 
massless spinor in the background of the two-charge black 
hole in AdSs; the Fermi momenta are kp = q — n — 1/2 
(in units of the chemical potential), where q is the charge 
of the spinor, and n labels the Fermi surface [19 . The 
RN-AdS black hole background is more complicated due 
to the extra feature of the oscillatory region. 

It is helpful to understand the normal modes by look- 
ing at the poles of the Green's function at arbitrary u. 
The numerical calculations suggest the following features, 
as illustrated in Fig. |4] At q — and fc = 0, all the poles 
of the Green's function are in the lower half cj-plane, and 
not close to = 0. As we increase the charge q, there 
are more and more poles moving across the origin to the 
upper half oj-plane, and the first one is labeled by n = 0. 
Suppose we start from a large q. As we increase fc, the 
poles will move to the right, across the origin to the lower 
half w-plane. Therefore, the largest fc corresponds to the 
ri = mode, which triggers the onset of the instability. 
The situation is similar to the fermionic case with a cru- 
cial difference that the quasibound states in the upper 
half plane in Fig. |4] are now in the lower half plane [20] . 

According to Ref. [B] , we can associate the first type of 
instability with the hybridized critical point, and the sec- 
ond type of instability with the bifurcating critical point. 
Recall that the system has rotational symmetry. In the 
fermionic case, a nonzero kp gives a spherical Fermi sur- 
face. However, in the bosonic case, the condensation of 
the scalar cannot happen at a spherical shell of the mo- 
mentum space. Instead, the quantum phase transition 
happens at w = and fc = 0. Without tuning extra pa- 
rameters, only the second type of instability can happen. 



5 



m =-4 



k 



2 3 
1 



k 




k 



FIG. 2. Phase diagram for the standard quantization. The sohd hnes correspond to zero modes. The osciUatory region is 
shaded, and will move to the right as we increase rn? . In the right plot, the tip of the oscillatory region corresponds to a 
bifurcating critical point at A: = 0. The dotted line is the BPS bound for q < A+/-\/3, where q is the R-charge [171 118j . 
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FIG. 3. Phase diagram for the alternative quantization. The oscillatory region is shaded. The dashed lines are solutions to 
Eq. (201, but do not represent zero modes. The dotted line is the BPS bound for q < A-/\/S. 



TABLE I. The critical value of the superfiuid velocity (± is for the standard and alternative quantization, respectively) 





q 


Superfiuid velocity 




-4 < < -3^ 


0<q< ^(1± Vm2 + 4) 
q > ^(1± +4) 


IR geometry instability at 
Zero mode instability at Sx ~ 


= kmJq 
ks/q 




0<g< 


No instability 




> -3 


^"''2+3 <q< ^/(l + Vm2+4) 


IR geometry instability at Sx 


= km/q 




q > ^(1 + +4) 


Zero mode instability at Sx ~ 


ks/q 



For the alternative quantization, the = —4 and —3 cases arc not included. 
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FIG. 4. Schematic plot of the poles of the Green's function 
when q is large. 



By tuning an extra parameter, both types of instability 
can happen at A; = 0. 



B. Bifurcating critical point 

When the IR scaling exponent i/^ at fc = becomes 
imaginary, there is a bifurcating critical point 6 . We 
will denote i' := i'k=o = ^^'sA*' where u = m? — 2q^ +3. 
From u > to u < 0, the Green's function at = 
keeps finite, but bifurcates into the complex plane and 
has a cusp [6]. 

Near = 0, the quantities Oj^-* and in the Green's 



6 



function can be expanded as 



,(0) 



r,(0) 



By Eq. (19 1, we have 



a ± i/feo; - 
^+^^r(A+-i) 



(28) 



V^(^^r(^ + 7s)r( 
r(A+-i)(^(- 



V2(m2+4)3-^+^+/'r(^ 



-^)+27^1nl8) 



V3J 



^r(A- 



-^j+27-lnl8j 



^2(m2+4)3" 



l + A_/2f 



(29) 



where i^{x) :— r'(a;)/r(a;) is the digamma function, and 
7 is the Euler-Mascheroni constant. It can be checked 
that aP - I3a = -l/{2\jTin? + 4) is satisfied. 

Near the bifurcating critical point j/ — > 0, the IR 
Green's function can be written as 



(30) 



where Qo{ui) is the IR Green's function when = 0, which 
can be obtained by Eq. (|8| with = 0: 



g^{uj) = - ln(-2iw) - 27 - 

Particularly, for a neutral scalar, 

( iuj 

go(w) = -ln(^-y 



V6 



7- 



(31) 



(32) 



By expanding the Green's function, Eq. (151, at small 



we obtain the Green's function for the bifurcating critical 
point: 



G 



/3go(c^)+/3 



(33) 



The finite temperature generalization of Q^i^uS) is given 
by Eqs. (|B9]) and (|B10[ ) in Appendix [b] 

Note that Eq. (331 is obtained by taking the v Q 
limit at fixed w. If we want to examine the poles of the 
Green's function, we need a Green's function valid to 
arbitrarily small uj. To do so, we need to use tJfc=o('^) 
instead of Qq{ijj). In the condensed side (u < 0), there 
are infinite number of poles in the upper half w-plane; 
these massive states will condense and lead to instability 

i- 



where Sx is the superfluid velocity, which is the source, 
and is the supercurrent, which is the expectation value 
PIH^ . The Klein-Gordon equation for a scalar cj) cou- 
pled with both At and sd lo — and fc = is 



m 
T2" 



Z2/ 



= 0. (35) 



When the superfluid velocity is too large, the superfluid 
phase will return to the normal phase. At the critical 
point, Ax = 5^; is a constant. From Eq. (35), we can 



see that the superfluid velocity plays the same role as 
the momentum fcj^ which is consistent with the fact that 
the superfluid velocity is the gradient of the phase of the 
order parameter. The phase diagrams are Figs. |2] and [3] 
by replacing k with qSx- 

According to previous studies of holographic supercon- 
ductors, there are two mechanisms for instabilities: (1) 
When the charge of the scalar is large, the effective mass 
mlfi = + g*-^q^A^ is negative enough to produce an 
unstable mode [H |28j. (2) When it is close to the zero 
temperature, and the effective mass m^g = [w? — 2q^)/12 
is below the AdS2 BF bound m|p = —1/4, the IR ge- 
ometry is unstable [3 13 US] ■ Here the quantum critical 
point can be reached by tuning the superfluid velocity, 
and the above two mechanisms for instabilities corre- 
spond to the two types of instabilities we discussed be- 
fore: zero mode instability and IR geometry instability, 
respectively. They can happen at the same time, giving 
a mixed critical point, as shown in Figs.[2]and|3j when a 
solid line touches the oscillatory region. 

We call the stable region {k or qSx > max(fcs, A;ir)) 
the normal phase, and we call the unstable region the 
superfluid phase. When the mass is small (—4 < m? < 
—3), there are two cases: (i) For small charge, the scalar 
is unstable due to the second type of instability; (ii) For 
large charge, the scalar is unstable due to the first type 
of instability. When the mass is large (m^ > _3)^ there 
are three cases: (i) For small charge, the scalar is stable; 
(ii) For intermediate charge, the scalar is unstable due 
to second type of instability; (iii) For large charge, the 
scalar is unstable due to the first type of instability. The 
result is summarized in Table HI 



C. Critical points with superfluid velocity 

A superfluid with a supercurrent flow can be studied 
holographically by turning on a vector potential 

Ax^Sx + + • • • , (34) 



^ We reinterpret k as the superfluid velocity (in tlie x-direction) . 
Another way to reinterpret k is the magnetic field added by Ay = 
Bx, assuming that the metric is not changed. After a separation 
of variables <I> ~ X{x)<p{z), the magnetic filed plays the same 
role as k'^ [qB <-> fc^) in the equation for ip l5l 1261 . In terms of 
a dyonic black hole in AdS4, the quantum critical point by the 
second type of instability is studied in Ref. |27| . 
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FIG. 5. A slight difference between AdS4 and AdSs. We set 
the AdS radius L = 1. In AdSs, we can only start from the 
standard quantization. 



IV. ADDING A DOUBLE TRACE 
DEFORMATION 

A. Hybridized critical point 

Without the superfluid velocity, the zero mode at A: = 
can also be achieved by tuning another parameter k^, 
which describes a double trace deformation in the bound- 
ary CFT: 



2 



(36) 



where (O) — B [6l|30]. The Green's function becomes 



(37) 



Similarly, starting from the alternative quantization, we 
can add a double trace deformation with coefhcient k_, 
which is related to k+ by k_ = The alternative 

quantization is allowed only if the of the scalar is in 
the interval < < fn^p + 1 [Uj, which leads to 
a slight difference between AdS4 and AdSs as illustrated 
in Fig. [Sj Since we are interested in the instability trig- 
gered by a zero mode at = 0, we do not want other 
instabilities to exist. Therefore, we require u > 0, i.e., 
the parameters are not in the oscillatory region, in which 
the system is already unstable. In AdS4, if we start from 
the alternative quantization, there is still an interval of 
not in the oscillatory region. In AdSs, however, we 
can only start from the standard quantizationj^ 
The Green's function at the leading order in uj is 



.{0} 



,(0) 



+ «+6(°)+gfe(L^)(a^"+K+&^'0 



(0) 



(38) 



The boundary condition for a pole in the Green's function 



at = is a 



(0) 



= 0. 



Thus, by Eq. (IT9|, the 



® When > —3, Eq. ( |36| is an irrelevant term. We assume that 
the UV geometry is not changed, and examine (both UV and 
IR) instabilities by the Green's function. 



critical value of k+ is 



A_-l 
2 



q 

V3 



ir( 



A_-l 
2 



- ^)r(A+ - 1) 



g 

V3 



ir(^ 



V3^ 



r(A_ - 1) 



(39) 



When = Kc, we obtain a hybridized critical 

point, which is described by an order parameter in the 
Ginzburg-Landau sector hybridized with a strongly cou- 
pled sector, the emergent CFTi dual to the IR AdS2 [6]. 
The Green's function near a hybridized critical point can 
be written as 



G{uj,k) 



1 



where 



(40) 



(0) 



C 



Vm2 +4(6^°))2 

iq_ 



r(2-)r(i-.-^) 



(41) 



For a neutral scalar, the first order in w vanishes, so we 
need to write h^uP' instead. 

We assume < 1/2 for a charged scalar, or j/ < 1 for 
a neutral scalar. For example, v < \ for a neutral scalar 
means —3 < < 9, which is already a large range. 
Then we do not need since the uP'^ term is dominant 
as a; — ?> 0. For the Green's function, w = is a branch 
point in the complex cj-plane, and we define the physical 
sheet as € (— 7r/2, 37r/2). The Green's function has a 
pole at 



hCiy) 

Note that /i < 0, and for q > 0, 



< arg 



1 



< 



V < 



(42) 



(43) 



(see Appendix [B]); for g = 0, C is real and negative. 
Therefore, when k+ < Kc, the pole will be in the upper 
half oj-plane of the physical sheet. When k+ > Kc, the 
pole will be either in the lower half w-plane of the physical 
sheet, or on a non-physical sheet. 

We can easily plot the critical value of k+ for the hy- 
bridized critical point as a function of u, as shown in 
Fig. [6j The boundary for the bifurcating critical point 
is w = 0. The two curves u = and — inter- 
sect at (m, K-I-) — (0,-2), which is a marginal critical 
point. When changes across the curve from k+ > Kc 
to < Kc, a pole will move across the origin to the 
upper half w-plane, causing an instability. We can see 
that Kc is a single- valued function of u, which implies the 
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following. If we keep increasing or decreasing after a 
pole moves across the origin of the complex w-plane, it 
will never come back to the origin. However, we cannot 
exclude the possibility that a pole can appear from infin- 
ity in the upper half w-plane, or a pole in the upper half 
cj-plane can move to infinity and disappear. The analytic 
solution cannot capture these features. We hope further 
numerical calculations can be helpful to make this clear[^ 
We will refer to the IR instability as the singularity 
of the Green's function at — 0, and UV instability 
as the singularity of the Green's function at \io\ — > oo 
in the upper half complex w-plane. For a CFT at zero 



density, the boundary for the UV instability is 







[5]. If the UV and IR instabilities are related, a tentative 
phase diagram for a neutral scalar is shown in Fig. |6] 
The boundary for the hybridized critical point is much 
more intricate than previous numerical results showed 
[B]- Especially, the curve will cross the u axis (k+ = 0). 
From Eq. (39), we can see that when 1 — A_ — — n, i.e., 
4, where n 



= (n +1) 



0, 1, 



the denomi- 



nator of Kc is infinity; thus, if there is no other infinity 
in the numerator, we have Kc = 0, which implies that 
the static Green's function G{uj = 0) — > oo. How ever , 



when A+ (or A_) is a integer, we need to use Eq. (A5) 



in Appendix [X] to calculate the Green's function, and the 
result is usually finite. This puzzle is largely related to 
the renormalization near integer values of Aj- (see Ap- 
pendix Ia|) . 



B. Marginal critical point 

The marginal critical point happens at both m = and 
K+ = Kc- Recall that the IR scaling dimension 5^ = 
1/2 -I- Vk\ when v — the double trace deformation is 
marginally relevant in the IR CFT 6] . The only possible 
parameters for the marginal critical point are m? = —3 
and q — Q. The critical value of is Kc = —2, which can 
be obtained from Eq. (39) by setting (7 = and fc = 
first, and then taking the — > —3 limit. The direct 
calculations are as follows. The solution of (p in the outer 
region at w = and fc = is 



Ci + In 



2z2 + 1 



The solution in the inner region is 



ln[12(l - z)] 
^12(1 - z) 



^12(1-^)' 



(44) 



(45) 



^ There are two challenges in numerical calculations. One is that 
the result becomes inaccurate very quickly as the Im(aj) becomes 
large; the other is that it is more difficult to obtain the expecta- 
tion value accurately when the scaling dimension is large. 



where Qq is given by Eq. ( 32 ) . By matching 
we obtain 



and 



1 



In 18 



2^6 
1 

2^6 



In 18 



/6 



(46) 



By expanding Eq. (38 1 at small the Green's function 



near the marginal critical point is 

/3eo(w)+/3 



G = 



{a + K+l3)go{uj) 



K+/3 



which gives 
G 



lncj-21n6-|-6 4-7-i7r/2 



(k+ 2)(lna; - 21n6 + 6 -^7 - i7r/2) - 12 



(47) 



(48) 



where the critical value of k+ is —2. 

We can replace the zero temperature Green's func- 
tion tjfc(a;) with the finite temperature Green's function, 
Eq. ( |B7[ ) in Appendix iBj and the phase diagram of T-k+ 
can be obtained [64. At finite temperature T << n, the 
result is 



G = 



(/^++2)[^(i-2^)+lnf 
whose imaginary part is 



T 
18 



imG ~ — tanh — - = 
24 2T 



TTU! 

48r 

TT 

.24 



sgn(a;) 



6-^7] - 12 

uj«T 
uj » T. 



, (49) 



(50) 



V. DISCUSSION 

Starting from the extremal RN-AdSs black hole, we 
have studied the quantum critical points by solving the 
Klein-Gordon equation in the bulk. The result gives us 
a glimpse of some strongly interacting systems, whose 
properties are beyond the reach of the perturbative 
method in quantum field theory. Just like the harmonic 
oscillator and the hydrogen atom as exactly solvable 
models capture essential features in quantum mechanics, 
the exactly solvable model in this work, together with 
a previous fermionic one [19], captures many essential 
features in AdS/CMT. 

We have calculated the Green's function to the leading 
order in uj. The non-analyticity of the Green's function 
indicates two types of instabilities: one is triggered by a 
zero mode, and the other is triggered by the instability 
of the IR geometry. In the standard/alternative quanti- 
zation, the zero modes of the system are always at finite 
k. However, we can tune an extra parameter to make 
the zero mode be at k = 0. We considered the quantum 
critical points of two systems, whose finite temperature 
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— »- bifurcatin 
marginal 


g \ 


- 


\— hybridized 





2 4 6 8 

U 

FIG. 6. The tentative phase diagram for a neutral scalar [q - 
consistent with Ref. [S]. The bifurcating critical point is at w = 
point is at k+ = k^. The curve crosses the horizontal axis at u 
and the blue shaded region when k+ > is for UV instability. 



K+ 




10 20 30 40 50 

u 

0), with different ranges oi u = rn? + 3. The small u part is 
0, and the oscillatory region is < 0. The hybridized critical 
- 8. The red shaded region when k+ < is for IR instability, 
'he white regions are stable. 




normal 




Sec. Ill C 



Sees. Ill B and IV 



FIG. 7. The left plot is the phase diagram for Sec. [mCl The 
parameter Sx is the superfluid velocity. SC denotes (holo- 
graphic) superconductor. The right plot is the phase diagram 
for Sees. |III B and IV according to Ref. [6]. The parameter g 



is u for the bifurcating critical point, or k+ for the hybridized 
critical point. 



phase diagrams are illustrated in Fig. [7| The extra pa- 
rameter in the first system is the superfluid velocity. The 
RN-AdSs geometry describes the normal phase, and de- 
pending on the parameters (m^,(?), the system can de- 
velop zero mode instability, IR geometry instability, or 
be stable. The extra parameter in the second system is 
the double trace deformation. The zero mode instability 
gives a hybridized critical point. In the second system, 
besides the above IR instabilities, there is also UV insta- 
bility, which is not captured by the analytic solution. 

There are several remaining questions in the second 
system as follows, (i) What is the full phase diagram in- 
cluding both UV and IR instabilities? We already have 
the boundary for the IR instabilities. To find the stable 
region, we need not only the boundary for the UV insta- 
bilities, but how the poles move as we change the param- 
eters, (ii) What happens to the Green's function when 
A+ approaches an integer? This is related to the curve 
crossing the u-axis (/t+ = 0) in Fig. |6] It seems that the 
A_(_ approaching an integer limit is not the same as the 
result obtained by setting A+ be the integer, (iii) What 
are the analogous models in condensed matter physics? 



Note added: Recently there is another related work 
IST for the two-charge black hole in AdSs, in which the 
massless Klein-Gordon equation is analytically solved. 
We can obtain that the zero modes as = q — 2n ~ 
where vi^ = ^/k'^ + \ and n is a nonnegative integer. Only 
the first type of instability exists, because the electric 
field approaches to zero in the near horizon limit |19j . 
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Appendix A: Mathematical notes 

Hypergeometric function 2-^1(0?, /3; 7; x). We will de- 
note 2F1 hy F for simplicity. The following formulas are 
helpful to obtain and understand the analytic solutions: 



Fia, /?; 7; x) = {l- x)-^F{-, - a, f3; 7; -) . (Al) 

X — I 

If Re(7) > Re(a-t-/3), 

r(7)r(7 - a - /?) 



F(a,/3;7;1) 



(A2) 



r(7 - a)r(7 - /?) 
A connection formula for the hypergeometric function is 

F{a, /3; 7; x) = AF{a, l3;a + p - -f + l;l - x) 

+ B{1 - xy~"~'^F{-f - a, 7 - /3; 7 - a - /? + 1; 1 - x), 

(A3) 
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where 



r(7)r(7 - a - /9) 

r(7-a)r(7-/3)' 



B = 



r(7)r(a + /3 - 7) 



r(a)r(/3) 



(A4) 



Other solutions for the Klein-Gordon equation ([7|. If 
A_|_ is an integer, Eq. ( 16 ) is no longer a general solution. 
If 2vk is not an integer, the general solution for Eq. ([T]) 

'At u ~ Q can be written as 



(z) =Ci 



Z^+(l~z2)-l/2-''. ^/A+-l 



(2z2 + l)-l/2-i.fc+A+/2 



2^1 



2 

A+-1 



V3' 2 



g 1 — ' 



9 



' (2z2 + l)-l/2+'^-c+A+/2 - 'V 2 V3' 2 

If 2g/%/3 is not an integer, the general solution can also be written as 

1 g A+ - 1 



1 q ^ 1- 

— h t^t H = ; 1 + 2i/i, ; — 

V3 '"2z2 + l 



2^^-A+/2+q/V3 .^^ 
2Fl\ 



(2z2 + 1)9/73 



[l-z') 



2->-A+/2-g/\/3 



(2^2 + l)-g/v^ 



A+-1 



V3' 
73' 



1 



2g 2z^ + 1 

J_ , , 2q_ 2z^ + 1 
^/3' ^/3' 1-^2 



(A5) 



(A6) 



r 



If all A+, 2vk, and 2q/^/i arc integers, it is likely that the 
solution is an elementary function. It is especially conve- 
nient to use Eq. ( |A5[ ), because the boundary condition at 
the horizon to obtain zero modes is simply Ci = 0. We 
will discuss some special cases, in which we only consider 
the standard quantization. 



Special case 1: w? = —A (A+ = 2). By using Eq. ( |A5[ ), 
the result for the static Green's fmiction is 



G{^ = Q) = ^[^[-^ + Vk 



Vk 



— ) 
V3J 



27 + In 3 , (A7) 



where ipix) is the digamma function. 



Special case 2: 



-3 (A+ — 3). Similarly, the 



result for the static Green's function is 
Giu = 0) = (3^.,^ - g2) [v.(l + - -|) 



J^k 



— ) 
V3J 



Special case 3: rn^ — (A+ 
solution for Eq. ([7]) at w = is 

<j)o = Ci+G2 I — ^ + ln 



4) and q = Q. The 



1 



1-Z2 



2z2 + l 



(A9) 



The solution for Eq. Q in the inner region is 0/ = 
gj^/[i2(i-2:)] obtain G(w = 0) = 0, which implies 

Kc — > oo at M = 3, and thus is consistent with Fig. |6] 

Special case 4- rn"^ = 5 (A+ = 5). This corresponds to 
u = 8 in Fig. [6] The result is cumbersome, but it is clear 
that the Green's function at w = is a finite number. 



which implies that Kc ^ 0. However, Fig. [6] shows that 
Kc = at tt = 8. This inconsistency can be explained 
as follows. If we take the limit that A+ approaches an 



integer, Eq. ( 16 1 diverges, and we need a renormalization. 



which changes Kc- 

Let's take a more careful examination about the cases 
when A+ = 2, 3, • • • . When A+ = 2, i.e., = m|p = 
—4, (f) — Az^ In z + Bz^ + • • • , where A is the source. 
Now we consider A+ > 2. According to Eq. (19), the 



expectation value B is divergent when A+ is an integer. 
If we add a small number 8 to A+, we have 



B = -+S, + 0(<5), 





(AlO) 



where B^ is the rcnormalized value of B. By the expan- 
sion 



z'' = l-H<51nz+i((51nz)2-h. 



(5^0, (All) 



27 + ln3 - ij + - - 3i^fe. (A8) we have 



Az^-(1 
Az^-{\ 



Bz^+^\\ + ---) 



+ + & In z + 0(<5 In z)) (1 + ■ 



(A12) 



where "• • • " denotes higher-order terms in z. The rcnor- 
malized Green's function by the standard quantization 
is Gr — Br/ A. Since there is a logarithm, the near 
boundary expansion should be evaluated at a UV cut- 
off z = Luv {L\]\IL « 1). Therefore, there are two 
noncommuting limits: (5—5-0 and Luv — ^ 0. The Green's 
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function obtained by Eq. ( 19 1 makes sense only if the The IR Green's function at zero temperature is 



following condition is satisfied: 



6 -In 



« 1. 



(A13) 



Appendix B: IR Green's function 

The IR Green's function Gk{^) ioi the scalar field is 
given in Ref. 9,. We will briefly review the result first. 
The IR geometry described by 



/ = 12(1 -2)2, At=2V6{l-z). 



Define the AdS2 coordinate 



(Bl) 



12(1 - z) 



(B2) 



The solution to the Klein-Gordon equation with the in- 
falling boundary condition at the horizon is 



W^_^,,(-2*^C), 



(B3) 



where Wx^f_i{x) is a Whittaker function with the following 
asymptotic behavior: 



T4^A,p(x)^e-^/2x^(l + •••), |x| 



By expanding Eq. (B3| at a;(^ — > 0, we obtain 



(B4) 



(B5) 



r(-2i^fe)r(i + - ^) 
r(2r.,)r(i-^.fe-^) ^ 



(B6) 



The finite temperature generalization at T << /x (chem- 
ical potential) is 



tT 



iq_ 



Near the bifurcating critical point — ^ 0, 



(B7) 



(B8) 



where 



iq 



.2 76 2tiT) 
-27-ln(47rr)-?/;(^ 
Particularly, for a neutral scalar, 



V6 



(B9) 



t?f'(a.) = -V.( 



1 



ILjJ 

2^1 



) — 7 — ln(7rr). 



(BIO) 



For q > and v < 1/2, Eq. (43) can be proved by 
1 / T{2u)T{^ - ly - iq. 



(-C ) \ r(-2i^)r(i + - iq, 

= arg[cos 7r(i/ — iq,)] 

= arctan(tan7ri/tanh7r(7*) < tti>. 



(Bll) 



where we have used the refiection formula r(a;)r(l — a;) = 
tt/ sinTTx, and the fact that we can put any positive num- 
ber inside arg(a;) without changing its value. 
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